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Abstract 

We show that the eigenvalues and eigenfunctions of the stargen- 
value equation can be completely expressed in terms of the correspond- 
ing eigenvalue problem for the quantum Hamiltonian. Our method 
00 '. makes use of a Weyl-type representation of the star-product and of 

the properties of the cross- Wigner transform, which appears as an 
intertwining operator. 
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1 Introduction and Motivation 



One of the key equations in the deformation quantization theory of Bayen 
et al. [H E] is, no doubt, the stargenvalue (for short *-genvalue) equation 
H * # = E*f> where * is the Moyal-Groenewold "star-product" [II [21 [7] . In 
this Letter we show that the *-genvalue equation can be completely solved 
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in terms of the usual eigenvalue/eigenfunction problem Hip = Eip where H 
is the Weyl operator with symbol H (and vice versa). The underlying idea 
is simple: we first rewrite the equation H -kfy = E^> in the form 

H(x + \ihd p ,p- \ihd x )^{x 1 p) = E^(x,p), 

where H(x + \ihd P) p — \ihd x ) is the Weyl operator with symbol 

U{z,0 = H[x-\C P ,p + \Q. 

We next show that the solutions of this equation and those of Hip = Eip 
can be obtained from each other using a family of intertwining operators 
(which is countable when H is essentially self-adjoin); these operators are 
up to a normalization factor, the cross- Wigner transforms ip i — > W(ip, cp) 
where (p describes the set of eigenfunctions of H . Our approach is inspired 
by previous work [4J of one of us on the time- dependent Torres- Vega [8] 
Schrodinger equation in phase space. 

Notation 

We will write z = (x,p) where ieR" and p G (M n )*. Operators 5(M n ) — > 
S'(R n ) are usually denoted by A,B, ... while operators S{R 2n ) — ► S'(R 2n ) 
are denoted by A, B, ... The Greek letters ip, (p, ... stand for functions defined 
on M. n while their capitalized counterparts $, ... denote functions defined 
on M 2n . We will make use of the symplectic Fourier transform defined for 
* G S(R 2n ) by the formula 

Ufa) = FpBiz) = (^) n / e-i<*^(z')dz> 

where a(z, z') = p-x' — p' -x is the standard symplectic form on W 1 x (W 1 )* = 
M. 2n (the dot ■ stands for the duality bracket; in practice p ■ x can be seen 
as the usual Euclidean scalar product under the identification (IR n )* = M n ). 
The symplectic Fourier transform is involutive: the identity on 

S'(R 2n ). 

2 Stargenvalue Equation: Short Review 

In view of Schwartz's kernel theorem every linear continuous operator A : 
<S(M n ) — > iS'(lR n ) can be represented, for ip G iS(R n ), in the form Aip(x) = 
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(JC A (x, •), V) with K A e S'(R n x R n ). By definition the contravariant (Weyl) 
symbol of A is the tempered distribution A defined by the Fourier transform 

a(x,p) = (er&U, Krfx + |(-), x - §(•)) . (1) 

Assume that B : SiW 1 ) — ► S{R n ); then the product C = AoB exists and 
its Weyl symbol is given by the Moyal product 

a * b(z) = (^) 2n / 1 e^ u ' v) a(z + \u)b{z - \v)dudv. (2) 

J JM. n xR n 

The main observation we will exploit in this paper is the following: if we 
write a = H and b — ^ then we can write 

H*V = HV, (3) 

where 

H = H(x + \ihd p ,p — \ihd x ) 

is a certain pseudodifferential operator on iS(M 2n ) we are going to identify. 
Let us view the linear operator H : ^ i — > H * ^ on 5(lR 2n ) as a Weyl 
operator. Using formula (j2J), the kernel of H is the distribution 

Ksi*>v) = (^) 2n / ^ {u > z - y) H{z-\u)du, (4) 

hence using ([1]) and the Fourier inversion formula the contravariant symbol 
of H is 

M(z, 0=[ e^JC 5 (z + ± V ,z- ± V )d V . 

JR 2n 

Using (j4j) and performing the change of variables u = 2z + rj — z' we get 
K s {z + lv,z- \rj) = (^n) 2n e 2 ^ [ ei^H{\z')dz 

JR 2n 

setting H(^z') = H\ji{z!^ the integral is {2nh) n times the symplectic Fourier 
transform F^H\/2{—rj) = (Hy2) a (—rj) so that 

={^eT [ e^e^(H l/2 ) a (- V )d V 

JM. 2n 
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where J = I I is the standard symplectic matrix. Since the second 



—I y 

equality is the inverse symplectic Fourier transform of (Hi/ 2 ) a calculated at 
the point z + J(. We finally get 

M{z,0=H{x-l( p ,p+lQ (5) 

where we are viewing ( = (( x ,(p) as the dual variable of z = (x,p); this 
justifies formula ([3]) viewing if as the quantized Hamiltonian obtained from 
EI by the quantum rule 

(x,p) i — > (x + ~ihd p ,p— \ihd x ). (6) 



3 i7-Calculus 

There is another very fruitful way of interpreting the Weyl operators H = 
H*ty. Let us return to the expression (J2J) with a = H and b — ^; performing 
the changes of variable u = 2(z' — z) and v = zq this formula can be rewritten 

as 



e -H z °> z ')H{z')dz' 



e -H z ^(z-lz )dz . 



Observing that the integral between brackets is (27ch) n times the symplectic 
Fourier transform of H we can write this formula in the form 

= (£hT / H^z )f(z Mz)dz (7) 

where T(z ) is the operator defined by 

f{z^{z)=e-^ z ^{z-\z Q ). (8) 
Formula (jTJ) is strongly reminiscent of the representation 

W = (^eT [ H*(z )T(z )Wza (9) 

of a Weyl operator H in terms of its covariant symbol H% = F^H and the 
Heisenberg-Weyl operator 

T(z )ip{x) = e^ p °- x -^°- Xo) ij(x - x Q ), 
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except that T(z ) is allowed to act on functions of z and not only of x. This 
feeling is amplified when one notes after a straightforward calculation that 
the operators T(zq) obey the relations 

f(zo + zi) = e-^ Z0 ^f (z )T ( Zl ) (10) 
f( Zl )f(z ) = e-^ Z0 ^f(z )f( Zl ) (11) 

which are similar to those satisfied by the Heisenberg-Weyl operators. These 
facts suggest that T(zo,t) = e^T(zo) defines a unitary representation of the 
Heisenberg group. Let us prove this is indeed the case. For this we will need 
the linear mapping : S(R n ) — ► S(R 2n ) defined by 

W^> = (2irh) n/2 W(ilj, 0) (12) 

where denotes an arbitrary function in <S(IR n ) such that | \<f>\ \ L 2 = 1. W(ip, 4>) 
is the cross- Wigner distribution; we thus have explicitly 

WM*) = (^)" /2 / e-^{x + \y)<t>{x-\y)dy. (13) 
Jr™ 

In view of Moyal's identity 

the operator extends into an isometry of L 2 (W l ) onto a subspace H<p of 
L 2 (R. 2n ); we are going to see in a moment is closed in L 2 (R 2n ), but let us 
first give a formula for the adjoint of W^. We have 

W;*(z) = (^) n/2 f e^-y^(2y - x)^(y,p)dpdy (14) 

jR n 

(it follows from a straightforward calculation using the identity (W^I^l^r 2 ") = 

(^|w;*) L2(R n)). 

Proposition 1 The range Ti^ of W</, is closed, and hence a Hilbert space. 

Proof. Set P^ = W^W^ where W£ is the adjoint of W^; we have = P£ 
and Pcf,P£ = P^ hence P^ is an orthogonal projection. Since W^W^ is the 
identity on L 2 {R n ) the range of is L 2 (M. n ) and that of P^ is therefore 
precisely Ti^. Since the range of a projection is closed, so is Tt^. ■ 
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This result, together with formula fill I) shows that T(zq) and T(zq) are 
unitarily equivalent representations of the Heisenberg group H ra ; the irre- 
ducibility of the representation T(zq) : H n — > follows from von Neu- 
mann's uniqueness theorem for the projective representations of the CCR. 

Let us return to the operator H = H*. A straightforward calculation 
showing that satisfies the intertwining relations 

x * W^ip = (x + lihd p )W^ip = W^xtp) 

p * W^ip = (p - ±ihd x )W^ = W^-iW^) 

an educated guess is then that more generally: 

Proposition 2 (i) The operator intertwines the operators T(zq) and 
f(z ): 

W*(f (zo)VO = T(z )W^; (15) 

(ii) We also have 

HW lP = W^H and W^H = HW£. (16) 

Proof. Making the change of variable y — y' + Xq in the definition (|T3l) of 

Wtf, we get 

W^f{zo)M){z) = e~i< z ^W^{z - \z ) 
which is precisely ( fi~5l) . Applying to both sides of (jSJ), we get 

W+m = (^)™ / H*{z Q )W^[T(zMdz Q . 

JR 2 " 

and hence 

WtpHip = (^) n [ H*(z )[T(z )W^]dz , 

JR 2 ™ 

which is the first equality ([TBI) in view of formula ([7j). To prove the second 
equality it suffices to apply the first to W£H = (H*W^)*. ■ 

4 Spectral Results 

We will need the following result, which is quite interesting by itself: 
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Lemma 3 Let (<pj)j be an arbitrary orthonormal basis of L 2 (M. n ). Then the 
vectors = W^<pk form an orthonormal basis of L 2 (R 2n ). 

Proof. Since the are isometries the vectors form an orthonormal 
system. It is sufficient to show that if \I/ G L 2 (R 2ri ) is orthogonal to the 
family ($j,k)j,k (and hence to all the spaces H.^) then \l/ = 0. Assume that 
(^|3>jjfc)z,2(R2n) = for all j, fc. Since we have 

it follows that W^,^> = for all j since (<pj)j is a basis; using the anti- linearity 
of in we have in fact WT^ = for all <p G L 2 (R n ). Let us show that 
this implies that * = 0. In view of formula (fill) for the adjoint of the 
operator W£ has kernel 

**(v,p) = (&) n/2 e**^ ) <l>Q*y-*)- 

Let us fix x; the property = for all (p is then equivalent to $ x ) = 
for all $ x G <S(R 2n ) (fixed 2) and hence * = 0, which we set out to show. ■ 
We now have everything we need to prove the main results of this Letter. 
We begin by stating the following general property: 

Theorem 4 The following properties are true: (i) The eigenvalues of the 
operators H and H = H* are the same; (ii) Let ip be an eigenf unction of H : 
Hip = Xip. Then, for every (p, the function ^ = W^ip is an eigenfunction of 
H corresponding to the same eigenvalue: H^f = A*. (Hi) Conversely, if^ is 
an eigenfunction of H thenip = Wl^S? is an eigenfunction of H corresponding 
to the same eigenvalue. 

Proof. That every eigenvalue of H also is an eigenvalue of H is clear: if 
Hip = Xip for some ip 7^ then 

H(W^iP) = W^HiP = A(W^) 

and W^ip 7^ because is injective; this proves at the same time that W^ip 
is an eigenfunction of H. Assume conversely that if \l/ = A\l/ for ^ 7^ and 
A G R. For every (p we have, using the second equality (fIBj) . 

Ml* = W1HV = XWl^f 
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hence A is an eigenvalue of H; is an an eigenfunction if it is different 

from zero. Let us prove this is indeed the case. We have W$WT§! = P^ 
where P^ is the orthogonal projection on the range H.^ of W$. Assume that 
= 0; then P^l = for every (f> G S(M n ), and hence ^ = in view of 
Lemma [3] above. ■ 

Remark 5 The result above is indeed quite general, because we do not make 
any assumption on the multiplicity of the (star) eigenvalues, nor do we assume 
that H is essentially self-adjoint. Notice that the proof actually works for 
arbitrary <fi G iS'(!R ra ). We present some examples at the end of this section. 

Corollary 6 Suppose that H is an essentially self-adjoint operator on L 2 (R n ) 
and that each of the eigenvalues X , Ai, Xj, ... has multiplicity one. Let 
i/)o,ipi, ... be a corresponding sequence of orthonormal eigenf unctions. 

Let tyj be an eigenfunction of H corresponding to the eigenvalue Xj . Then 
there exists a sequence (atj jk ) k of complex numbers such that 

*i = Yl a i^it mth *i>* = W +$i e H i n H t- ( 17 ) 

i 

Proof. We know from Theorem H] above that H and H have same eigenvalues 
and that ^j )k = W^ k ipj satisfies the eigenvalue equation H^/j^ = Xj^j )k . 
Since H is self-adjoint and its eigenvalues are distinct, its eigenfunctions ipj 
form an orthonormal basis of L 2 (R n ); it follows from Lemma [3] that the tyj jk 
form an orthonormal basis of L 2 (R 2n ), hence there exist non-zero scalars ctj^,i 
such that tyj = i a j,k,e^k,e- We have, by linearity and using the fact that 

HV kjt = Xk9 kjt , 

On the other hand we also have H^fj = Xjtyj, 

which is only possible if atj,k,i = for k ^ j; setting a^i = ctjj/ formula ffTTl) 
follows. (That tyj t t G TCj D Tie is clear using the definition of Tig and the 
sesquilinearity of the cross- Wigner transform.) ■ 
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We remark that the continuous spectrum can be dealt with in a simi- 
lar fashion provided that one generalizes the transform by allowing the 
"parameter" to be a tempered distribution (in which case the normalization 
condition ||</>||l2 = 1 does no longer make sense, of course); the same remark 
applies to the case where H is no longer essentially self-adjoint (cf. the re- 
mark following the proof of Theorem H]). To illustrate this, let us consider 
the two following typical examples (in dimension n = 1): 

• H(x,p) = p. In this case H = —ihd/dx is a symmetric operator and 
the equation Hif) = Eip has solutions for every real value of E; these 
solutions are the tempered distributions ip(x) = C exp(iEx/h) (C any 
complex constant). A straightforward calculation shows that 

W^{x,p) = C'e% {E - p)x F0(pj 

where C is a new constant and Fcf) is the Fourier transform of <fi. If we 
let <p range over <S'(M n ) and use the fact that the Fourier transform is 
an automorphism of <S'(M n ) we see that W^ip can be any distribution 
of the type 

*(x,p) = §(p)e^ {E - p)x 

with $ e iS'(lR n ); these distributions are precisely the solutions of the 
stargenvalue equation 

p * \J> = (p - ±ihd x )$! = Ety 

clS Si straightforward calculation shows. 

• H(x,p) = x. Here H is the operator of multiplication by x; this a 
symmetric operator without any eigenvalues and eigenf unctions. It is 
however self-adjoint, and the solutions of Hip = Eip are the distribu- 
tions ip = C5(x — E); one finds by an argument similar to that above 
that W^ip can be any distribution of the type 

#(x,p) = $(x)e-T(^-*)f 

which is the general solution of the stargenvalue equation 

x * $ = (x + \ihd p )^ = E^. 
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The previous treatment of the stargenvlue equation for operators with a 
continuous spectrum can be made rigorous in the setting of Gelfand triples. 
In our setting (S(IR n ), L 2 (R n ), S'(R n )) is the Gelfand triple of interest and 
the corresponding weak formulation of the eigenvalues and eigenvectors of an 
operator from <S'(K. n ) to iS(M n ). For further information on Gelfand triples 
we refer the reader to the standard reference [3J. 

5 An Example and its Extension 

As an illustration consider the harmonic oscillator Hamiltonian 

H= l -(p 2 + x 2 ). (18) 

In view of the results above the spectra of the operators H and H are iden- 
tical. Choosing for simplicity fi = 1 the eigenvalues of H are the numbers 
A at = N + ~ with N = 0,1,2,.... The normalized eigenf unctions are the 
rescaled Hermite functions 

Mx) = (2 k k\V^)^e-^ 2 n k (x). (19) 

where 

H k (x) = (-l) km e x2 (£) k e- x2 

is the k-th Hermite polynomial. Using definition (112j) of W$ together with 
known formulae for the cross-Wigner transform of Hermite functions (see for 
instance Wong [10J, Chapter 24, Theorem 24.1) one finds that the eigenfunc- 
tions of H are linear superpositions of the functions 

* j+k)k (z) = (-iy ( I -^) l 2f +1 C ?£ 4(2|z| 2 )e-^ 2 (20) 

where ( = x + ip and ^jj+ k = ^j+ k , k for k — 0,1, 2, here 

C){x) = jjx- k e x {±Y {e~ x xi +k ) ,x>0 

is the Laguerre polynomial of degree j and order k. (For similar results see 
Bayen et al. [2].) 

Notice that the example above can be generalized without difficulty to 
the case of arbitrary quadratic Hamiltonians of the type 

H = -Mz ■ z 
2 
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where M is a positive-definite symmetric matrix. In fact, in view of Williamson's 
diagonalization theorem there exists a symplectic matrix S such that 



M = S T DS , D 



A 
A 



where A is the diagonal matrix whose entries are the moduli Uj > of the 
eigenvalues ±iujj of JM. We thus have 

3=1 

and H o S = SHS^ 1 where S is anyone of the two metaplectic operators 
associated with S. The eigenvalues of H o 5 and are the same; they are 
the numbers 



and then the eigenfunctions ips of H o S and those, tp, of if by the formula 
ips — Sip. Now, the eigenfunctions of H are tensor products of rescaled 
Hermite functions; using the fact that ips = Sip together with the symplectic 
covariance formula 

W(SiP,S<P)(z) = W(iP,<P)(S- 1 z) 

satisfied by the cross- Wigner distributions, the eigenvalues of H = H-k are 
calculated in terms of tensor products of the functions (1201) . We do not 
give the details of the calculations here since they are rather lengthy but 
straightforward. 



6 Concluding Remarks 

Due to limitation of length there are several aspects of our approach we 
have not discussed in this Letter. For instance, he methods developed here 
should apply with a few modifications (but in a rather straightforward way) 
to more general phase space (for instance co-adjoint orbits). A perhaps 
even more exciting problem is the following, which is closely related to our 
previous results [5] on the relationship between the uncertainty principle 
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and the topological notion of symplectic capacity. A rather straightforward 
extension of the methods we used in [5J shows that if 

\W^(z) \ < Ce-^ (aN2+%|2) for z G M 2n (21) 

then we must have ab < 1. In particular we can have IW^VK' 2 )! < Ce - ^' 2 ' 2 
only if e > fi\ it follows that the Hilbert spaces Ti^ do not contain any 
nontrivial function with compact support: assume in fact that ^ G Ti^ is 
such that ^{z) = for \z\ > R > 0. Then, given an arbitrary e < h one 
can find a constant C £ such that 1^(2)1 < C £ e~~\ z \ 2 , which is impossible 
since \1/ = W^ip for some ip G L 2 (IR n ). This suggests (taking Theorem H] 
into account) that the solutions \1/ of the *-genvalue equation cannot be too 
concentrated around a point in phase-space. In fact we conjecture that if 
an estimate of the type \^{z) \ < Ce^^ Mz ' z (M symmetric positive-definite) 
holds for an eigenfunction of the stargenvalue equation, then the symplectic 
capacity of the ellipsoid Mz ■ z < h must be at least \h. We will come back 
to this topic in a near future. 

Acknowledgement 7 The authors would like to express their deep gratitude 
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